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Relativistic mean field calculations suggest that the charge and matter deformations significantly 
differ in some of the unstable neutron and proton rich nuclei. We discuss the effects of the difference 
on the fusion reactions induced by them at energies near and below the Coulomb barrier by taking 
the ' ' Na + Pb reactions as examples. We also discuss whether one can probe the difference 
by the so called fusion barrier distribution analysis. 
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£ '. I. INTRODUCTION 

Secondary beam experiments are opening up new physics to be studied for a variety of unstable nuclei across the 
\Q ' beta-stability line in the nuclear chart. One of the basic questions is to study their charge and matter deformations. 
We are especially interested in the difference between the charge and matter deformations in each unstable nucleus. 
One cannot, however, apply the standard techniques such as the gamma ray spectroscopy which have been used 
for stable nuclei. A possibility is to use heavy-ion fusion reactions at energies below or near the Coulomb barrier, 
which have recently been established to be very sensitive to the shape of the target and/or projectile nuclei |I]J^]. An 
alternative is to predict the charge and matter deformations theoretically based on the non-relativistic Hartree Fock 
or the relativistic mean field calculations with various refinements such as including the pairing correlations. 
\Q ' In this paper, we take Na isotopes as an example to anticipate the variation of the charge and matter deformations 
along an isotope chain and discuss whether the significant difference between them expected for some isotopes can 
be probed through the analysis of the fusion reactions induced by them by assuming 08 Pb as the target nucleus. 
We use the relativistic mean field calculations to assess the charge and matter deformations for each Na isotope. We 
then calculate the fusion cross section a(E) based on the coupled-channels formalism assuming thus obtained values 
of deformations and taking the ground state rotational band into account. We then calculate the associated so called 
q ' fusion barrier distribution d ^np^ §|. 

The paper is organized as follows. In sect. 2 we describe the formalism of our coupled-channels calculations, i.e. the 
simplified framework based on the no-Coriolis approximation, somewhat in detail, because the formalism is not well 
known for the case of fusion of a deformed odd mass projectile which we are going to discuss in this paper. In sect. 
3, we present the results of the relativistic mean field calculations and of the coupled-channels analyses. The paper is 
summarized in sect. 4. 
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II. SIMPLIFIED COUPLED-CHANNELS FORMALISM FOR THE FUSION OF A DEFORMED 

ODD-MASS NUCLEUS 

We discuss in the next section the fusion of deformed odd mass Na isotopes with 208 Pb target. We calculate the 
fusion cross section by taking the excitation of the ground state rotational band of the Na projectile into account with 
the coupled-channels formalism in the no-Coriolis approximation. In this section, we describe the major aspects of 
the method. 

We denote the coordinates of the relative motion between the projectile and target nuclei and of the rotational 
motion of the projectile, i.e. a Na isotope, by R and £, respectively. We ignore the intrinsic excitation of 208 Pb. The 
total Hamiltonian is then given by 
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H(R,t)=T R + Vo{R)+H (Q + V(R,£) (1) 

with 

f = - — — — R 2 — + J^—L 2 
R 2fiR 2 dR dR 2y,R 2 

Vq(R) is the sum of monopole nuclear and Coulomb interactions between the projectile and target nuclei. We call it 
the bare interaction. Hq(^) is the intrinsic Hamiltonian of each Na projectile, whose eigenfunctions and the eigenvalues 
are denoted as ff mi (£ ) and cik , 

ffo(OvL,(0 = e/*PjU(0 (2) 

We assume that each deformed Na isotope can be well described by the strong coupling model 0] for the rotational 
excitation, where the wave function is given by 



vLAZ) = \m&r{ Dl m I KfflxK + (-lJ'+^-xCtyx-*} (3) 

On the r.h.s. of eq. ([?]), we used ft for £ in order to explicitly indicate the physical meaning as the Euler angles, xk 
is the intrinsic wave function. We assume the deformation to be axially symmetric, and K is the component of the 
angular momentum of the rotational motion along the symmetry axis. 

The V(R, £) in cq.(Q) represents the interaction between the projectile and target nuclei after the monopole part has 
been subtracted. We assume the same form of coupling Hamiltonian as that is familiar in coupled-channels calculations 
for the fusion reactions of even-even nuclei. It consists of the nuclear and Coulomb parts, and is expressed as 

V(R,0 = f N (R)Y,Yx(R) -TMO+^/f (R)Y Xli (R)T^(0 (4) 

in terms of the spherical harmonics for the relative motion and the tensor operators 7\ M which cause the rotational 
excitation of the Na projectile. We further assume the latter to be given by 

Tx»(0 = £ OAvi>^(n) = OAoX^cn) = /3Ax£„(n) (5) 

v 

where a\ v are the deformation parameters in the body fixed frame Q. We have kept only the v — term to be 
consistent with the assumption of the axially symmetric deformation. The deformation parameter f3\ is identified 
with the nuclear and the charge deformation parameters, (3^ and in the nuclear and Coulomb coupling terms, 
respectively. The f N (R) and f£(R) are the nuclear and Coulomb coupling form factors, which we assume to be, 

jy_ , d — Vq p iZpZ^e 2 Rp 



r(R) = -Rp-m U , ,„, D D1 /„ ]■ fnR) 



dR L l + exp[R~ i? ]/ao 2A + 1 R^ 1 

with Rq = Rp + Rt, Rp and Rt being the nuclear radii of the projectile and target. 
As usual, we introduce the channel wave functions by 



iJM 
^IKL 



(R0= ^2 {LmLlrni\JM)YL mL (R)tpf mi {£) (6) 



and expand the total wave function for a fixed total angular momentum J and its z component M as, 

* JM (R,Z) = ]T ^j KL {R)^ L {R,0 (7) 



R 

IKL 



It is then straight-forward to obtain the following coupled 2nd order differential equations which determine the radial 
wave functions of the relative motion &'} KL (R) , 

Vo(R) + ^ m L(L + l) + e IK -EU-} KL (R)+ ]T Vikl,i>k>l>(R)H K 'l>(R) = (8) 



2fi dR 2 y ' 2(iR 2 



XI'K'L' 
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where the coupling matrix elements are given by 



V, 



IKLJ'K'L 



(R) 



f N (R)Px +fx(R)Px 



(-1) 



L+J+l 



l'+i'A + (21' + 1)(2L + 1)(2A + 1)(2L' + 1) 



4-7T 



L A V 




L I J 
V L' X 



(X0I'K\IK)S KK , (9) 
we used the explicit expression of the tensor operator given by eq.(pj), which 



In obtaining eq.(|8|) with eq.(| 
leads to the selection rule Skk 1 - 

The full coupled-channels equations (|^) become fairly intricate if many channels are included. In this paper, we 
use the no-Coriolis approximation, in other words the rotating frame approximation, which considerably reduces the 
dimension of the coupled equations ||. The process to calculate the fusion cross section then consists of three major 
steps. The first is to replace the channel dependent centrifugal potential energy in eq.(^) by that in the entrance 
channel, 



L{L + l)h 2 Li(Li + l)h 2 



2^R 2 2[iR 2 

where Li is the initial orbital angular momentum. 

The second step is to perform the unitary transformation, 



(10) 



which transforms the original coupled-channels equations eq.(||) into 

U{U + i)h 2 



2[i dR 2 



V (R) 



2fiR 2 



(ii) 



(12) 



with 



I+2X-I' 



f N {R)P\ + fxi R )Px 



.l + (-l) A 1 /( 2A + 1)(2J' + 1) 
4tt(2I+1) 



(XOI' K\IK)(X0I' - K\I - K) (13) 



The K is the projection of the intrinsic spin of the Na projectile along the z-axis of the rotating frame which is taken 
to be parallel to R. Eq.(|lJ) shows that only states with the same K quantum number couple to each other. The 
dimension of the coupled-channels equations is thus considerably reduced. 

The third step is to solve the reduced coupled-channels equations for each K=—Ii, —Ii + 1, — 1,7,, Jj being 
the ground state spin of the Na projectile. Note that K = Ii is fixed. As usually done for heavy-ion fusion reactions, 
we replace the proper boundary condition that the radial wave functions for the relative motion should be regular at 
the origin R = by the incoming wave boundary condition at an absorption radius R a bs inside the potential pocket, 
and require 



rphiK 

ik ex P 



^l b {k LdKk {R'))dR' 



Hj(ki K R)5i U + R L .f 'H+{k IK R), 



, R < Rabs 
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with 



ki K = j2[i/h\E -e IK ) 



iw 



2fiR 2 



V Q (R) 



(14) 

(15) 
(16) 



In eq.@, Tf-i K and R^ are the transmission and reflection coefficients, respectively. We used these notations 



IK 



instead of and Rj^ because all the transmission and the reflection coefficients are the same for a given Li 

irrespective of J. Once they are determined, the barrier transmission probability is given by 
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P jR (E) = p( J ) L * R (E) = khaKk 



k 



rpLiK 

IK 



(17) 



where k = kj^ is the wave number in the entrance channel. The first equality in eq.(|17|) means that the barrier 
transmission probability is the same for all J for a given Li. 

The fusion cross section can be related to p( J ) L i K by calculating the flux at R a b s - To that end, we use the inverse 
transform of eq. (11), 



^kl(R) = \j 2 ^\^ LQlk \ jk ) ¥ ik ( 18 ) 
k 



Averaging over the initial and summing over the final spin projections, we obtain 

7T 1 



^ = F2FTl£ E W+VP [J)LiK iE) (19) 

K,L z J=L,-h 

(2Lt + l)P L > k {E) (20) 



k 2 



K.L, 



We have used a simplified notation P LiK (E) in eq.(|2fj) because what one actually does is to solve the coupled equations 
eq. (|l^) for each set of K and Li and determine the barrier transmission probability. Eq. (po| ) reduces to the well known 
formula 



J 

for the fusion between two even-even nuclei by setting the initial intrinsic spin L to be zero. 



III. THE CHARGE AND MATTER DEFORMATIONS OF NA ISOTOPES AND THEIR EFFECTS ON 

THE FUSION REACTION 

We now discuss the charge and matter deformations of Na isotopes and the influence of their difference on fusion 
reactions taking 19 > 25 > 37 Na+ 208 Pb scattering as examples. We perform the relativistic mean field (rmf) calculations 
with the NLSH parameters j^J^] to learn the r.m.s. radii and deformation parameters of the neutron, proton, charge 
and matter distributions of Na isotopes. Though the pairing correlation may strongly influence the accurate values of 
deformation @], we ignore it, because our interest in this paper is to discuss global behaviors. The results are shown 
in Fig.l, where the solid circles, solid squares and solid triangles represent the r.m.s. radii of the proton, matter and 
charge distributions, respectively. The open diamonds are the effective root mean square matter radii deduced from 
the measured interaction cross section using a Glauber- type calculation || . We see that the theoretical results given 
by the solid squares well follow the general trend of the experimental data given by the open diamonds. Fig. 1(b) shows 
the isotope variation of the quadrupole deformation parameters of the proton, matter and neutron distributions. It 
indicates that the charge and matter deformations are significantly different for many isotopes far away from the beta 
stability line. 

We now discuss how this difference influences the fusion reactions, or conversely whether one can probe this difference 
through the analysis of the fusion cross section. We choose the 19 > 25 - 37 Na isotopes as representatives and analyze their 
fusion reactions with the 208 Pb target. As already mentioned, we treat 208 Pb as a structureless particle. 

Before discussing the results of the coupled-channels calculations, we show in Fig. 2 the potential barriers corre- 
sponding to fixed orientations of the deformed Na projectile assuming the quadrupole coupling, i.e. A = 2, 

V(R, 9) = ~ V ° + £^ + fC(R)P 2 (co S 9) (21) 

1 + exp[(R -R - R P l3?^P 2 (cos8))/a } V 47r 

where the potential parameters are assumed to be (Vo, ao) = (114.2MeV, 1.222/m). The radius parameters Rp and Rt 
have been calculated from the r.m.s. radius given by the rmf calculations according to Rp(Rt) — ^/| x r.m.s. radius 
of the projectile(target). The higher and lower barriers in Fig. 2 are for 70.12 and 30.55 degrees, respectively. These 
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are the effective angies in the so calied orientation average formalism in the degenerate spectrum limit approximation 
when one discusses the effects of the ground state rotational band on the fusion of a deformed even-even nucleus by 
truncating at the first excited 2 + state |l(J . The dotted line is the potential barrier when the effect of deformation 
is ignored. For each angle, the solid line is the potential barrier obtained by using different nuclear and Coulomb 
deformation parameters suggested by the relativistic mean field calculations. It is compared with the dashed line, 
which has been obtained by ignoring the difference between the charge and matter deformations. It is interesting 
to observe that the difference of the matter and charge deformations modify the barriers in opposite directions for 
proton and neutron rich isotopes. 

We now present the results of the coupled-channels calculations. There are no experimental data yet on the 
excited states of 19 > 25 > 37 Na isotopes. As we stated before, we assume a strong coupling rotational model, where the K 
quantum number is assigned to be K=| for all of them from the Nilsson diagram using the deformations suggested 
by the relativistic mean field calculations. The excitation energy is given by 

,,<= W + % KV (22, 

where J is the moment of inertia. 

It is known that the experimental moment of inertia is smaller than the rigid body value J r i g by about a factor 
of 2 . Though the deviation of the moment of inertia from the rigid body value is somewhat smaller for odd mass 
nuclei Q), we assume, 

J ~ Jrig/^ (23) 

for the moment of inertia of 19 < 25 > 37 Na nuclei. The rigid body value is given by 

J rig = ^AMR 2 (l + ^6) (24) 

where S is the deformation parameter, which is related to the quadrupole deformation parameter (32 as 

5 = 0.945/3 2 [l - ^ir 2 {a /Rp) 2 } + 0.34/3 2 . (25) 

In eq. (p5|), the correction due to the surface diffuseness given through the diffuseness parameter ao has been included 
only to the leading order. The quadrupole deformation parameter 02 is identified with the matter(nuclear) deformation 
parameter fi^ °f the Na isotopes. 

We solve the coupled-channels equations taking up to the I m ax=7/2 member of the ground state rotational band 
and including only quadrupole deformation for all three reactions. The fusion cross section is then obtained by eq. 
( f2p| ) . Once the fusion excitation function has been obtained, the fusion barrier distribution is calculated by the point 
difference formula of Ai?=2MeV in the center of mass energy. The results are shown in Fig. 3. Each figure contains 
three theoretical lines. The dotted line has been obtained by the potential model. The solid line has been obtained by 
the coupled channels calculations using the charge and matter deformation parameters obtained from the relativistic 
mean field calculations. The dashed line has been obtained in the same way, but by ignoring the difference between 
the charge and matter deformations. We used the value of (3 N for both nuclear and charge deformations in this case. 

These figures show that the effects of the difference between the charge and matter deformations are significant 
in the 37 Na+ 208 Pb fusion reactions. They are noticeable also in the 19 Na+ 208 Pb fusion reactions. The existence of 
three or more peaks can be expected because of the truncation at I max =7 /2 by considering the degenerate spectrum 
limit. There appear indeed three separated peaks for the 37 Na+ 208 Pb reaction. The shift of the highest and lowest 
peak positions by the difference between the charge and matter deformations accords with that shown in Fig. 2 for 
the two channel problem. Unfortunately, no visible effects are seen for the 25 Na+ 208 Pb fusion reactions, which would 
be much more tractable experimentally than the other reactions. Though 25 Na has a large difference in the charge 
and matter deformations comparable to that in 19 Na and 37 Na, its effect on the fusion cross section is compensated 
by the counter effect of the difference between the charge and matter radii. 

The above conclusions might depend on the special choice of the moment of inertia given by eq. ( |23| ) . In order to 
test this, we repeat the same calculations for 19,37 Na+ 208 Pb reactions by changing the moment of inertia to 

J « J rig /3 (26) 



In Fig. 4, the dashed and the solid lines are the same as those in Fig. 3 obtained by assuming eq. (|23|). These lines 
are compared with the short dashed and the dash-dotted lines which have been obtained by assuming eq. (^6|). The 
figures show that the fusion barrier distribution distinguishes the different choice of the moment of inertia and that 
it clearly reflects the difference of the charge and matter deformations for both choices of the moment of inertia. 
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IV. SUMMARY 



The isotope variation of the charge and matter deformations is an interesting question, which will be addressed 
with the advent of secondary beam experiments. In this paper, we used relativistic mean field calculations to suggest 
that the charge and matter deformations can significantly differ from each other for some nuclei far away from the 
beta stability line. We then performed coupled-channels calculations using these informations, and have shown that 
the fusion excitation function and the fusion barrier distributions noticeably reflect the difference between the matter 
and the charge deformations in 19,37 Na+ 208 Pb reactions. Unfortunately, this effect is compensated by the opposite 
effect of the difference between the matter and charge radii in the 25 Na+ 208 Pb fusion reactions which will be more 
easily tractable experimentally. 

Finally, we wish to make some comments on the limitations of our relativistic mean field calculations. We ignored 
the pairing correlation in our relativistic mean field calculations of Na isotopes. The bulk properties of the light 
nuclei through helium to oxygen has been calculated by the relativistic mean field theory with and without pairing 
correlation p| . The study shows that the r.m.s. radii of light nuclei are not so sensitive to pairing correlation, whereas 
the deformation parameters are affected considerably. We are now extending the relativistic mean field calculations 
to include the effects of pairing correlation in order to have a more reliable estimate of the deformation parameters. 
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Figure Captions 



FIG. 1 

The ground state properties of Na isotopes calculated by the RMF theory without pairing correlation. The top panel 
shows the proton, matter and charge r.m.s. radii. The open diamonds show the measured matter r.m.s. radius taken 
from || . The bottom panel shows the quadrupole deformation parameters of the proton, matter and neutron density 
distributions. 

FIG. 2 

The effective fusion barriers for the s-wave scattering of 19 ' 37 Na from 208 Pb as functions of the separation distance R 
between the projectile and target nuclei. The dotted line is the bare potential. The solid and dashed lines have been 
obtained by respecting the difference between the charge and matter deformations, and by ignoring it, respectively. 
For each system, the upper and the lower barriers are for the 70.12 and 30.55 degree orientations, respectively. 

FIG. 3 

Effects of the difference between the nuclear and Coulomb deformations on (a) the fusion excitation functions and (b) 
the fusion barrier distributions in the 19 > 25 > 37 Na+ 208 Pb reactions. The dotted lines are for the potential model. The 
solid lines represent the results when the difference of the nuclear (matter) and the Coulomb(charge) deformations 
is taken into account. The dashed lines have been calculated by ignoring the difference of the charge and matter 
deformations. 

FIG. 4 

The sensitivity of (a) the fusion excitation function, and (b) the fusion barrier distribution to the choice of the moment 
of inertia in the 19,37 Na+ 208 Pb reactions. The dashed and the solid lines are the same as those in Fig. 3. They are 
compared with the short dashed and the dash-dotted lines obtained by assuming a different value for the moment of 
inertia in the coupled-channels calculations. 
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Fig. 1 
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